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ABSTRACT 

We present Sloan Digital Sky Survey and United States Naval Observatory 
observations of the V1647 Ori protostar and surrounding field near NGC 2068. 
V1647 Ori, the likely driving source for HH 23, brightened significantly in Novem- 
ber 2003. Analysis of SDSS imaging acquired in November 1998 and February 
2002 during the quiescent state, recent USNO photometry, and published 2MASS 
and Gemini data shows that the color changes associated with brightening sug- 
gest an EXor outburst rather than a simple dust clearing event. 
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1. Introduction 

In January 2004, J.W. McNeil discovered a new reflection nebula in the dark cloud Lynds 
1630 near M 78 (McNeil 2004) This object, now known as McNeil's Nebula, is apparently 
associated with an EXor-type eruption (Reipurth & Aspin 2004) of the embedded protostar 
V1647 Ori. 

EXors belong to the class of pre- main- sequence optical variables (Herbig 1977). They 
are classical T Tau stars which undergo irregular outbursts in the optical/UV of several 
magnitudes, named for the prototype EX Lup (Herbig et al. 2001). These outbursts are in- 
terpreted as episodes of substantial mass transfer resulting from instabilities in the accretion 
disk; they are present very early in the evolution of a protostar, as shown by the detection 
of EXor outbursts from deeply embedded Class 1 protostars in the Serpens star formation 
region (Hodapp et al. 1996). 

Clark (1991) first identified V1647 Ori as the young stellar object IRAS 05436-0007 on 
the basis of its IRAS colors. I-band and [SII] narrow-band imaging of the region by Eisloffel 
& Mundt (1997) revealed a faint I band source at the position of the IRAS object and 
reflection nebulosity extending to the north, identifying V1647 Ori as the likely driver for 
HH 23, located 170 arcsec north of the star. The bolometric flux of the source derived from 
IRAS and sub-millimeter photometry by Lis, Menten, & Zylka (1999) yields a luminosity of 
2.7 Lq and an inferred molecular gas mass of 0.4 M© assuming a distance of 400 pc to the 
Orion star- forming complex (Anthony- Twarog 1982). 

Estimates of the extinction towards V1647 Ori, Ay = 11™ - 15™, are found to be similar 
from photometry taken during the quiescent phase (Abraham et al. 2004) and during the 
eruptive phase (Reipurth & Aspin 2004; Vacca et al. 2004; Briceno et al. 2004; Andrews, 
Rothberg & Simon 2004). Thus, it is not clear whether the appearance of McNeil's Nebula is 
due only to the eruption of VI 647 Ori or to the eruption plus additional clearing of obscuring 
circumstellar dust. In this paper, we examine this question using pre-eruption multiband 
optical and near infrared data from the Sloan Digital Sky Survey (SDSS) and the Two Micron 
All Sky Survey (2MASS) compared with post-eruption data in the SDSS and 2MASS bands 
observed at the United States Naval Observatory. 



2. Observations 

We detect the protostar at four epochs of Sloan Digital Sky Survey (SDSS) imaging 
as a point source (SDSS J054613. 14-000604.1) coincident with the 2MASS X-band position 
(Q;2ooo=05'^46™13.1^ 52ooo=-00°06'05"). The SDSS observations consist of two pairs of scans 
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acquired in November 1998 and February 2002. Figure 1 shows an SDSS composite image 
made with the riz filters in which the protostar and the faint nebulosity to the north can be 
seen. 

A technical summary of the SDSS is given by York et al. (2000). The SDSS imaging 
camera is described by Gunn et al. (1998). The Early Data Release and the Data Release One 
are described by Stoughton et al. (2002) and Abazajian et al. (2003). The former includes 
an extensive discussion of the data outputs and software. Pier et al. (2003) describe the 
astrometric calibration of the survey and the network of primary photometric standard stars 
is described by Smith et al. (2002). The photometric system itself is defined by Fukugita 
et al (1996), and the system which monitors the site photometricity by Hogg et al. (2001). 
Abazajian et al. (2003) discuss the differences between the native SDSS 2.5m ugriz system 
and the u'g'r'i'z' standard star system defined on the USNO 1.0 m (Smith et al. 2002). 

The SDSS low Galactic latitude data which includes the Orion equatorial imaging used 
in this work are described by Finkbeiner et al. (2004). 

The U.S. Naval Observatory Flagstafi^ Station 1.0 m and 1.55 m telescopes were used to 
obtain eruptive phase BVRIz' and JHK photometry of V1647 Ori. For BVRIz', frames 
were taken and fiatfielded using twilight sky flats. DAOPHOT PSF fltting as implemented 
in IRAF was used to obtain photometric measures of the target since there is a bright part 
of the nebula only a few arcsec distant. For each dataset, ensemble differential photometry 
was performed using a set of secondary standard stars calibrated on two photometric nights 
with the USNO 1.0 m telescope. The z' measures were relative to z' secondary standard 
stars calibrated by the SDSS PT telescope (Hogg et al. 2001). For JHK, we used 2MASS 
stars in the fleld, eliminating obvious variables, to calibrate the data. We used a standard 
K filter for 2004 February 11, but used a K' filter for 2004 April 12. 



3. Results 

The pre-eruptive and recent optical and near-infrared photometry of V1647 Ori are sum- 
marized in Tables 1,2, and 3. The quiescent and eruptive phase spectral energy distributions 
are shown as uFi, in Figure 2. 

We examine the evolution of reddening invariant colors from the quiescent to the eruptive 

states to constrain the underlying physical processes. These colors take the generic form 
Cxyz = (x — y) — (y — z) * E(x — y)/E(y — z) where x, y, and z are the observed magnitudes 
in each passband. A color change is computed as ACxyz, where color changes having ACxyz 
statistically distinct from indicate changes in the spectral energy distribution (SED) not 
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consistent with pure dust clearing. The conversion from E(B — V) to extinction in each band 
follows Schlegel, Finkbeiner, & Davis (1998) and Finkbeiner et al.(2004), in preparation 

The reddening invariant colors measured before and during the eruption are listed in 
Table 4 for both Ry — 3.1 and Ry = 5.5, the latter appropriate for the larger dust grains 
found in star formation regions. Ry, defined as Av/E(B — V), is the ratio of the general to 
selective extinction. 

Characteristic grain sizes are inferred to increase from 0.17 /im to 0.21 /im as Ry ranges 
from 3.1 to 5.5 assuming a mix of silicate and carbonaceous populations (Weingartner & 
Draine 2001). In the coldest portions of the molecular clouds additional grain species such 
as refractory and volatile organics, olivine, water ice, orthopyroxene, trolite, and metallic 
iron are expected to contribute to the opacity (Pollack et al. 1994). Vacca et al. (2004) 
detect absorption due to water ice at 3.1 /^m in the near-IR outburst spectrum. 

Study of the selective extinction in background stars has revealed that Ry is not con- 
stant within the interiors of dark molecular clouds. Strafella et al. (2001) find that for the 
Bok globule CB 107 Ry increases inwards reaching a value of 6.5 at the core. Given the 
complex nature of the V1647 Ori environment including outburst cavities, accretion disks, 
and remnant envelopes it is likely that Ry can vary on small spatial scales. 

3.1. Determination of Extinction 

In the foregoing analysis we choose to characterize the intervening dust in terms of the 
SDSS z band extinction and Ry. Since the shape of the extinction curve is dependent upon 
the grain size distribution for A < O.dfim (Cardelli, Clayton & Mathis 1989) the traditional 
Ay and E(B — V) quantities are sensitive to both the dust column density and Ry. By 
adopting reference passbands longward of 0.9 /xm we remove the dependency on the form of 
the dust law although Ry must still be considered in the optical and UV. Alternate passbands 
that are free of Ry effects include the near-IR J and Johnson-Cousins /, the latter used by 
Weingartner k Draine (2001) who adopted Ai/N(HI) = 2.5 x 10~^^ cm^ for a standard 
gas-to-dust ratio. 

Following Abraham et al. (2004) and Reipurth & Aspin (2004) we assume that the 
quiescent phase near-IR colors are that of an embedded low-mass Classical T Tauri and 

obtain E(J — H) = 1.43™^ by dereddcning onto the Classical T Tauri locus of Meyer, Calvct, 
& Hillenbrand (1997). This value of E(J-H) corresponds to Aj = 3.39™ (A, = 6.40™) given 
Aeff(J) = 1.25/Lim and Xet[{H) = 1.65/xm and application of the Cardelli, Clayton & Mathis 
(1989) methodology. 
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3.2. Changes in the Spectral Energy Distribution 

In order to analyze the general differences between the quiescent and eruptive appear- 
ance we combine sets of observations to form representative values. The quiescent state is 
taken as the October 7 (2MASS J,H,K) and the 1998 November 17 (SDSS r,i,z) obser- 
vations. The optical photometry for the eruptive state is obtained from the 2004 February 
14 Gemini r and i data (Reipurth & Aspin 2004), and the 2004 February 26 USNO z band 
observations. The J, H, and K data for the eruptive state are from the 2004 February 11 
USNO observations. 

Table 5 presents the reddening invariant colors in the quiescent and eruptive states for 
Ry = 3.1 and Ry = 5.5. As evidenced by the greater than 5a change in all colors except 
for Criz, which has a marginal r band detection in the quiescent state, it is clear that the 
increase in emission is not consistent with a dust clearing event and thus we conclude that 
an intrinsic change occurred to the SED of V1647 Ori. 

In Figure 3 we compare the observed magnitude variations against AA^ — —2'^ and 
_4"^ (Just clearing events and, as concluded above based on reddening invariant colors, find 
they are not explainable in terms of a simple diminishing of the line of sight extinction. 

3.3. The Quiescent Phase 

From observations of Class II protostars (Classical T Tauris) we would expect that the 
veiling continuum due to magnetospheric accretion will be present in the SDSS u and g 

bands but weakening into the redder r, i and z bands (Calvet & GuUbring 1998). In a 
similar fashion, the thermal emission from the inner portion of the circumstellar disk should 
be bright in the H and K bands but diminishing into the bluer J band (Meyer, Calvet, & 
Hillenbrand 1997). For all but the mostly heavily veiled stars the i and z bands will only 
be affected by extinction. None of the reddening invariant colors can be used to give the 
intrinsic spectral type. 

Young (1 Myr) T Tauri stars exhibit spectral types between MO and M4 for masses of 
0.1 to 1.0 Mq (Baraffe et al. 1998). The use of the Meyer, Calvet, & Hillenbrand (1997) 
near-IR Classical T Tauri locus by Abraham et al. (2004) and Reipurth & Aspin (2004) 
imphcitly assumes a similar spectral type as this locus is based on observations of K7/M0 
stars in the Taurus star formation region. The i — z colors for these early to mid M spectral 
types range from 0.38 to 0.80 with considerable scatter on the order of several tenths of 
magnitudes (West et al. 2004). Given the observed i — z (2.01±0.06), dereddening to these 
intrinsic colors requires A^ = 3.1"^ to 4.2"* for Ry = 3.1 and 4.3"* to 5.7"* for Ry = 5.5. 



-6- 



Dereddening using the estimated ~ 6.4 and Ry = 3.1 results in an intrinsic i — z of 
-0.47, which is too blue for the stellar locus. Assumption of Ry = 5.5 yields an i — 2; of 0.20, 
corresponding to a late K spectral type (Finlator et al. 2000). If V1647 Ori is indeed a low 
mass protostar with a late K or M spectral type than either Ry is high or significant i band 
veiling is present in the quiescent state. 

In Figure 4 wc compare the observed quiescent phase SED against an MO photosphere 
seen under an extinction of Az=6.4 magnitudes. Excess emission is seen in the J, H, and K 
bands that could indicate the presence of a circumstellar disk. 

3.4. The Eruptive Phase 

The optical spectra acquired during the outburst by Walter et al. (2004) lack the TiO 
molecular absorption bands characteristic of late K and M stars. Walter et al. (2004) at- 
tribute this to cither an early photosphcric spectral type or overwhelming veiling. We adopt 
the latter interpretation given the apparent change to the intrinsic SED and the large bright- 
ness increase in the optical (Figure 3). 

If the new component to the SED is due to an EXor outburst then we expect to see 
emission from the high temperature (6000 - 8000 K) inner disk (Bell et al. 1995). In Figure 5 
we show the observed flux (z/F,^) increase in comparison with a 7000 K blackbody reddened 
by Az = 3.2™ and 6.4"^ for Ry = 3.1 and 5.5. In common with Walter et al. (2004) we see 
that a single EXor-like high temperature component can not reproduce the observed excess 
in both the optical and the near-IR. 

For the partial dust clearing events suggested by Rcipurth & Aspin (2004) and Walter 
et al. (2004) we find that due to the increase in optical depth towards the blue the additional 
flux due to a new source must fall more rapidly at the shorter wavelengths. The observed 
change in J — if is -0.61"*, which if entirely due to dust clearing requires AA^ = -3.2"*. The 
intrinsic flux increase required in this case is shown in Figure 5 and would require either a 
lower temperature for the new source or a significantly higher extinction than anticipated. 

Spectral indices are commonly defined by the ratio of the flux of the feature (F^) and 
that of the nearby pseudo-continuum (Fc), so the intrinsic measure of the feature strength 
is lo — Fg/Fc (Gizis 1997). If a veihng continuum is present, defined by F^ — rFc, then the 
measured spectral index, 7, is 
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Figure 6 shows the observed depression of the hne relative to the continuum level (1 — /) 
on a logarithmic scale as a function of the veiling, r, for unveiled indices of Jq = 0.0, 0.1, to 
0.9. From Figure 3 we see that the flux in the i band has increased by Am = -5 or a factor 
of ~100 which means the TiO features will be visible as less than a 1% fluctuation in the 
observed continuum. 



The observed photometric variations are not consistent with a simple dust clearing event 
as evidenced by changes in reddening invariant indices. We infer that the process of eruption 
involves changes to both the optical and near-IR SED. 

Application of an inferred prc-outburst extinction of = 6.4™ suggests a late K spectral 
type if the r — i and i — z colors are minimally affected by the veiling continuum arising 
from magnetospheric accretion shocks and that Ry has a value expected for a star formation 
region. An early to mid M spectral type is possible if the i band includes an excess emission 
of 0.2"* to 0.6"*, corresponding to a veihng between 0.2 and 0.7. 

We interpret the quiescent phase of V1647 Ori as an embedded low mass Classical 
T Tauri. Comparison with a reddened MO photosphere shows a near-IR excess that is 
presumably due to thermal emission from a normal circumstellar disk. 

The outburst SED is dominated by the new component. We see that a single blackbody 
having the ~7000 K temperature expected for an EXor inner disk can not simultaneously 
reproduce both the optical and near-IR portions of the SED. It is possible, as suggested 
by Reipurth & Aspin (2004), that a partial dust clearing event occurred in combination 
with an intrinsic brightening. We note that a reduction of the line of sight extinction would 
preferentially lessen the required flux increase at shorter wavelengths. This steepening of 
the curve shown in Figure 5 would increase the difficulty of fitting the high temperature (B 
spectral type) component observed by Briceiio et al. (2004). 

Further study of V1647 Ori planned using the facilities at USNO, Apache Point Ob- 
servatory, and the Spitzer Space Telescope may clarify the nature and evolutionary state of 
this young star. 



4. 



Conclusions 
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Table 1. riz Photometry of V1647 Ori 



Date 


r 


i 


z 


Telescope 


1998 Nov 17 


23.04 ±0.22 


20.81 ±0.05 


18.80 ±0.04 


SDSS 


1998 Nov 28 


24.00 ± 0.80 


21.03 ±0.09 


19.19 ±0.08 


SDSS 


2002 Feb 07 


22.69 ±0.29 


20.77 ±0.08 


18.77 ±0.06 


SDSS 


2002 Feb 09 


23.88 ± 0.73 


21.33 ±0.11 


19.33 ±0.07 


SDSS 


2004 Feb 14 


17.4 


15.6 




Gemini^ 


2004 Feb 26 






14.60 ±0.06 


USNO 1.55m 


2004 Apr 16 






14.79 ± 0.06 


USNO 1.55m 


2004 Apr 27 






14.72 ± 0.09 


USNO 1.55m 



^From Reipurth & Aspin (2004). 
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Table 2. Johnson- Cousins Photometry of V1647 Ori 



Date 


B 


V 


Rc 


Ic 


Telescope 


2004 Feb 11 




18.52 ±0.07 


16.91 ±0.08 


14.92 ±0.05 


USNO 1.0m 


2004 Feb 16 


20.39 ±0.12 


18.59 ±0.06 


16.87 ±0.07 


14.90 ± 0.04 


USNO 1.0m 


2004 Feb 25 




18.68 ± 0.04 






USNO 1.55m 


2004 Feb 26 


20.84 ±0.08 


18.99 ± 0.04 




15.15 ±0.02 


USNO 1.55m 


2004 Apr 13 


20.74 ±0.33 


19.12 ±0.08 


17.24 ± 0.09 


15.28 ±0.05 


USNO 1.0m 


2004 Apr 16 


21.13 ±0.07 


19.26 ±0.03 




15.37 ±0.01 


USNO 1.55m 
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Table 3. Near-IR Photometry of V1647 Ori 



Date 


J 


H 


K 


Telescope 


1998 Oct 07 


14.74 ± 0.03 


12.16 ±0.03 


10.27 ±0.02 


2MASS 


2004 Feb 03 


11.1 ±0.1 


9.0 ±0.1 


7.4 ±0.1 


Gemini'^ 


2004 Feb 11 


10.79 ±0.01 


8.83 ±0.01 


7.72 ±0.01^ 


USNO 1.55m 


2004 Apr 12 


10.94 ±0.02 


9.06 ±0.02 


7.59 ±0.02*^ 


USNO 1.55m 



'"From Reipurth & Aspin (2004). 

•^The USNO 1.55m observations used standard K on 2004 February 11 and K' on 2004 
April 12, both placed on the 2MASS zeropoint system, but probably with transformation 
differences. 
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Table 4. Reddening Invariant Colors 



Color Rv = 3.1 Rv = 5.5 



Criz (r-i)- 0.987 * {i - z) {r - i) - 1.002 * {i - z) 

Q^j (i-z)- 0.824 * {z-J) (i-z)- 0.605 * {z - J) 

C,jH (z-J)- 2.443 *{J-H) (z-J)- 2.443 * {J - H) 

CjHK {J - H)- 1.563 *{H-K) {J - H) - 1.563 * {H - K) 
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Table 5. Reddening Invariant Color Changes 



Color 


Rv 


= 3.1 




Rv 


= 5.5 




Quiescenf^ 


Eruptive*^ 




Quiescent^ 


Eruptive'^ 




Criz 


0.25±0.23 


0.81±0.10 


2.2 


0.22±0.23 


0.80±0.11 


2.3 


CizJ 


-1.34±0.08 ■ 


■2.14±0.09 


-6.7 


-0.45±0.07 ■ 


■1.31±0.09 


-7.7 


CzJH 


-2.24±0.12 - 


-0.98±0.07 


9.4 


-2.24±0.12 - 


-0.98±0.07 


9.4 


CjHK 


-0.37±0.07 


0.23±0.03 


8.0 


-0.37±0.07 


0.23±0.03 


8.0 



'"The quiescent state is defined by the 1998 October 7 2MASS and 1998 November 17 
SDSS data. 

''The eruptive phase data are from the 2004 February 14 Gemini r and i, the 2004 February 
26 USNO z, and the 2004 February 11 USNO J, and K measurements. 

'^No" is the ratio of color change to the quadrature sum of measurement errors. 
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Fig. 1. — SDSS pre-eruption riz band mosaic image of V1647 Ori. The location of 
V1647 Ori is marked on this 2x2 binned mosaic which maps the SDSS r, i, and z bands onto 
blue, green, and red. Hcrbig-Haro objects, such as the large HH 24 complex {bottom center), 
are seen in blue due to the Ha emission appearing in the r band. The image is roughly 10 
arcminutes on a side and is displayed using a negated grayscale. North is up and east is to 
the left. The faint emission immediately north of V1647 Ori significantly brightens during 
the eruptive phase where it is seen as McNeil's Nebula. The diffuse r band emission due to 
HH 23, which may be driven by V1647 Ori, is highlighted by the box in the upper center. 
A color rendition of this image is available in the electronic version of the journal. 
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Fig. 2. — Quiescent and Eruptive Spectral Energy Distributions. The optical and 
near-IR SEDs are shown for the quiescent state (squares) and during eruption (circles). Data 
acquired by the SDSS and at the USNO are indicated by filled symbols. The open squares 
arc the 2MASS JHK observations and the open circles are the r and i band eruptive phase 
measurements of Reipurth & Aspin (2004). 
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Fig. 3. — Changes in Spectral Energy Distributions. The observed magnitude differ- 
ences (squares) are compared against pure dust clearing events for AA^ = -2"* (circles) and 
AAz = -4"* (triangles) and for Ry = 3.1 (solid) or Ry = 5.5 (open) dust laws. 
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|Lim 

Fig. 4. — Modeling the Quiescent Spectral Energy Distribution. Observed pie- 
eruption SED compared with an MO photosphere (3800 K) seen under a line of sight red- 
dening of Az = 6.4"* for Rv = 3.1 {solid line) and Ry = 5.5 {dotted line). Excess emission 
is evident longward of 1.2 /im (J). 
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Fig. 5. — Modeling the Eruptive Component. This figure compares the observed 
flux increase {filled squares) against a 7000 K blackbody viewed under (Rv=3.1 (solid), 
5.5 (dotted)) extinctions of = 3.2"* and 6.4"*, scaled to match the observed change in 
the z band flux. The triangles show the required intrinsic flux increase in the case of a 
AAz = —3.2™ partial dust clearing for Ry — 3.1 (filled) and 5.5 (open). 
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Fig. 6. — Effect of Veiling on Spectral Features. This figure plots the observed depres- 
sion of an absorption line relative to the continuum level (1 — /) on a logarithmic scale as 
a function of the veiling, r, for unveiled indices of /q = 0.0, 0.1, to 0.9. The dashed lines 
indicate the eruptive phase i band veiling of 100 and a relative depression level of 1%. 
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